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Section- 1 

INTRODUCTION 

The  transition  metal  dichalcogenide  (MX2)  materials,  and  in  particular  MoSo,  have 
been  one  of  the  most  widely  studied  classes  of  solid  lubricants.*"^  In  general,  the  effective 
lubricating  properties  of  M0S2  have  been  attributed  to  the  highly  anisotropic  structure  of  this 
material  (Fig.  1-1). 


van  der  Waals 
interaction 


FIGURE  1-1.  Schematic  side-view  of  the  M0S2  structure  illustrating  covalently  bonded 
S-M-S  layers  that  interact  primarily  through  dispersion  forces. 

The  weak  interlayer  (S-Mo-S)/(S-Mo-S)  interactions  enable  facile  interlayer  shear,  and  hence 
lead  to  a  small  coefficient  of  friction.  However,  factors  in  addition  to  the  layered  structural 
motif  must  also  be  important  since  many  structurally  similar  MX2  solids  (e.g.,  NbSe2)  are 
poor  lubricants.  Other  factors  believed  to  play  a  role  in  determining  the  overall  effectiveness 
of  a  lubricant  include;  (1)  the  electronic  structure  of  the  solid;  (2)  adhesion  of  the  lubricant  to 
the  substrate  (metallic)  interface,  and  (3)  reactivity  and  wear  of  the  lubricant  surface.  In  the 
case  of  M0S2  and  NbSe2,  photoelectron  spectroscopy  data  and  molecular  orbital  analyses 
have  been  used  to  argue  that  electronic  structure  can  affect  significantly  the  interlayer 
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shear. Specifically,  long-range  interlayer  interactions  between  the  half-filled  Nb-dz- 
orbitals  in  NbSeo  may  inhibit  interlayer  shear  relative  to  M0S2  Mtematively,  it  is  also 
possible  that  the  electronic  properties  of  NbSe2  promote  deleterious  surface  reactions  (and 
thus  wear)  or  reduce  the  binding  interaction  to  metal  substrate  surfaces  compared  with  M0S2. 

To  develop  a  microscopic  understanding  of  MoSt  and  other  metal  dichalcogenides, 
which  is  essential  for  the  rational  design  and  development  of  new  materials  with  improved 
tribological  properties  we  have  used  scanning  tunneling  microscopy  (STM)  and  atomic  force 
microscopy  (AFM).^'^^  STM  and  AFM  provide  direct  atomic  resolution  data  addressing  the 
structural,  electronic,  wear  and  frictional  propenies  of  interfaces  and  as  such  these  new 
techniques  are  essential  for  tribology  research.  In  this  report  we  summarize  the  first 
systematic  application  of  these  methodologies  to  the  transition  metal  dichalcogenide 
materials. 
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Section-2 


STRUCTURE  AND  ELECTRONIC  PROPERTIES 
OF  METAL  DICHALCOGENIDES 

2.1  EXPERIMENTAL  METHODS 

The  measurements  made  during  this  project  were  carried  out  on  single  crystal  metal 
dichalocgenide  samples.  Although  single  crystals  are  experimentally  difficult  to  prepare,  they 
represent  the  best  defined  samples  from  which  strong  conclusions  can  be  drawn.  Facilities  for 
single  crystal  growth  of  metal  dichalcogenide  solids  were  developed  during  the  initial  stages 
of  the  project. 

Single  crystals  of  M0S2,  Ni0.1Mo0.9S2.  MoS2.xSex  (x  =  0.25,  0.5)  and  MoS1.75Teo.25 
were  grown  from  polycrystalline  powders  by  chemical  vapor  transport.22  For  example,  a 
stoichiometric  mixture  (5  g  total  mass)  of  molybdenum  and  sulfur  was  sealed  in  a  quartz  tube 
under  vacuum  and  reacted  for  10  days  at  1000°C.  Following  this  initial  reaction,  the  resulting 
polycr>  :talline  powder  was  ground  resealed  under  vacuum  in  a  quartz  tube  together  with 
excess  sulfur  (60  mg)  and  iodine  (100  mg).  Single  crystals  were  grown  in  a  50  -  100  °C 
gradient  over  a  three-week  period.  The  crystals  obtained  from  the  growth  region  had  the 
expected  plate-like  layered  morphology  and  varied  in  size  from  1  mm  x  1  mm  to  3  mm  x  3 
mm.  The  metal  (Ni)  and  chalcogenide  (Se,  Te)  doped  M0S2  crystals  were  grown  using 
similar  conditions.  The  procedures  for  the  growth  of  other  metal  dichalcogenide  samples 
examined  during  this  project  (e.g.,  NbSe2  and  TaS2)  were  also  similar  and  have  been 
published. 

Both  commercial  and  custom  STMs  have  been  utilized  in  our  measurements.  The 
underlying  basis  for  STM  experiments  and  its  application  to  metal  dichalcogenide  materials 
have  been  discussed  previously  in  several  of  our  published  reviews,  The  home-built 
ultrahigh  vacuum  (UHV)  STM  system  was  constructed  during  the  project  to  carry  out 
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experiments  under  highly  controlled  vacuum  conditions.  This  instrument  is  illustrated  in 
Figure  2-1. 


I-v.ipor.UDfs 


FIGURE  2-1.  Schematic  view  of  the  UHV-STM  system. 

Preliminary  studies  of  wear  and  reactivity  were  carried  out  with  this  instrument  and  are 
discussed  below.  In  addition,  we  anticipate  that  this  specialized  instrument  will  be  utilized 
extensively  in  our  new  AFOSR  project. 

In  addition,  we  have  developed  an  AFM  for  operation  in  controlled  environments.  In 
force  microscopy,  an  interface  is  imaged  by  scanning  a  surface  below  with  a  probe  that  is 
attached  to  a  cantilever;  deflections  of  the  cantilever  are  related  to  surface  features.  The 
force(s)  that  cause  the  cantilever  deflections  define  the  mode  of  imaging.  For  example,  when 
the  tip  is  in  contact  with  the  surface  repulsive  electrostatic  forces  usually  are  the  dominant 
interaction;  this  contact  regime  is  termed  repulsive  mode  imaging.  Alternatively,  if  the  tip  is 
removed  from  the  surface  a  van  der  Waals  attractive  interaction  may  be  the  dominant  force 
between  the  sample  and  tip;  this  case  is  termed  attractive  mode  imaging.  All  of  our  AFM 
studies  reviewed  below  were  carried  out  in  the  repulsive  mode  using  a  modified  commercial 
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instrument.‘'^-56  5131^4  cantilever-tip  assemblies  were  used  for  both  imaging  and  modification 
studies.  An  illustration  of  the  experimental  apparatus  is  shown  in  Figure  2-2. 


XYZ  Translator^ 

FIGURE  2-2.  Schematic  view  of  the  atomic  force  microscope. 

2.2  Ni-DOPED  xMoS2 

Our  initial  studies  of  the  atomic  level  structure  and  electronic  properties  of  solid 
lubricants  focused  on  M0S2  and  Ni-substituted  M0S2  NixMoi.xS2.^^  The  Ni-doped  materials 
are  a  particularly  important  system  for  high-resolution  studies  since  recent  work  has  shown 
that  Ni-substitution  enhances  the  tribological  properties  of  M0S2  compared  to  the  pure 
compound.  To  understand  these  results  we  have  used  STM  and  AFM  to  characterize  the 
structural  properties  of  well-defined  single  crystals  of  NixMoi.xS2.  First,  stoichiometric 
crystals  with  x(Ni)  =  0,  0.02,  0.05,  and  0.10  were  synthesized.  The  Ni-doped  materials  have 
been  characterized  by  a  variety  of  techniques  including  variable-temperature  resistivity 
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measurements  and  X-ray  diffraction.  These  studies  show  that  Ni-doping  does  not 
significantly  alter  the  average  structural  propenies  of  the  M0S2  system. 

AFM  images,  which  reflect  primarily  surface  structure,  exhibit  a  regular  atomic 
lattice  with  a  period  of  3.15  ±  0.05  A.  This  structure  corresponds  to  the  surface  sulfur  atom 
positions  and  is  in  agreement  with  the  structure  inferred  from  diffraction  measurements.  In 
contrast,  STM  images  of  the  Ni-doped  materials  exhibit  pronounced  defects  in  the  hexagonal 
atomic  lattice  (Fig.  2-3).  Since  the  AFM  data  show  conclusively  that  there  are  no  large 
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FIGURE  2-3.  30  x  30  gray-scale  STM  images  of  (a)  Nio.1Moo.9S2.  and  (b)  M0S2. 
Localized  effects  are  clearly  visible  in  (b). 


11 


structural  variations  due  to  Ni-substitution  these  defects  can  be  attributed  to  localized 
variations  in  the  density  of  electronic  states  due  to  nickel  ate  ms  in  the  lattice.  Hence,  we  can 
conclude  that  on  an  atomic  scale  Ni-substitution  only  affects  electronic  states  near  the  Fenui 
level  and  does  not  perturb  significantly  the  structure.  Such  localized  changes  are  not  e.xpected 
to  enhance  layer-layer  sliding  or  material  plasticity,  although  these  defects  could  enhance 
adhesion. 

2.3  ANION-DOPED  MoS: 

We  have  also  characterized  the  structural  and  electronic  effects  of  selenium  and 
tellurium  substitution  in  .VI0S2.  The  initial  motivation  for  this  work  was  the  suggestion  that 
Se  or  Te  dopants  could  function  as  atomic  scale  bearings  as  shown  schematically  in  Figure  2- 
4. 


Mo  Mo  Mo 


Tej  S 

Mo  Mo 


FIGURE  2-4.  Schematic  cross-section  of  M0S2  substituted  with  Se  and  Te.  The  larger  radii 
of  these  anions  is  highlighted. 

Alternatively,  it  is  expected  that  the  selenium  or  tellurium  dopants  might  function  as 
good  substrate  anchoring  sites. 

STM  and  AEM  investigations  of  these  new  materials  have  produced  several 
interesting  results.-^  Atomic  resolution  STM  and  AFM  images  of  MoSi,75Seo.25  exhibit 
similar  hexagonal  lattices  with  periods  of  0.32  ±  0.01  nm  and  0.32  ±  0.02  nm,  respectively 
(Fig.  2-5). 
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FIGURE  2-5.  9  nm  x  9  nm  (a)  STM  and  (b)  AFM  images  of  MoSj  vsSeoj.s-  The  white 
scale  bar  corresponds  to  0.9  nm. 

In  addition,  the  observed  atomic  corrugations  in  these  STM  and  AFM  images  are 
similar  to  the  corrugations  determined  from  images  of  M0S2.  These  STM  and  AFM 
experiments  have  thus  shown  that  the  structural  and  electronic  properties  of  M0S2  are  not 
perturbed  significantly  by  Se-substitution.  It  is  expected,  therefore,  that  the  coefficients  of 
friction  for  MoS1.75Seo.25  and  M0S2  will  be  similar. 

In  contrast,  we  have  discovered  significant  atomic  level  structural  changes  at  the 
surfaces  of  Te-substirated  M0S2.  Representative  STM  and  AFM  images  are  shown  in  Figure 
2-6. 
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FIGURE  2-6  8  nm  X  S  nm  (a)  STM  and  (b)  AFM  images  of  MoSi_75Teo.25.  The  small 
arrow  in  (b)  highlights  one  local  atomic  protrusion. 


The  STM  images  exhibit  hexagonal  structure  with  a  period  and  vertical  corrugation  similar  to 
that  determined  for  pure  M0S2.  These  results  have  shown  that  there  are  no  large  perturbations 
in  the  electronic  states  upon  Te-substitution.  The  AFM  images  have  shown,  however,  well- 
defined  atomic  size  protrusions  (Fig.  2-6  (b)).  Since  we  have  shown  that  the  number  of 
protrusions  .scales  directly  with  the  Te  concentration  it  is  likely  that  they  correspond  directly 
to  the  tellurium  atoms  in  the  lattice.  Hence,  it  is  possible  that  Te-substitution  may  reduce  the 
layer-layer  shear  energy  and  corresponding  sliding  friction. 
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Section-3 


REACTIVITY  AND  WEAR  OF  M0S2  AND  NbSe2 

3.1  WEAR  OF  M0S2  VERSUS  NbSe2 

The  effective  lubricating  properties  of  quasi-two-dimensional  materials  such  as  M0S2 
is  typically  ascribed  to  facile  interlayer  shear  of  the  van  der  Waals  bonded  layers.-^-'^ 

However,  factors  other  than  this  layered  structural  motif  must  also  be  important  since 
structurally  similar  transition  metal  dichalcogenide  materials  are  not  all  good  lubricants.  In 
the  case  of  M0S2  and  NbSe2,  the  inferior  properties  of  NbSe2  have  been  attributed  to  the  less 
facile  interlayer  shear  that  results  from  the  different  electronic  configuration  of  Nb(IV)  versus 
Mo(IV).‘^-5  Direct  support  of  this  proposal  and  more  generally  a  detailed  microscopic 
understanding  of  friction  and  wear  in  these  materials  have,  however,  not  been  available. 

Using  AFM  we  have  been  able  to  provide  new  and  general  data  addressing  the  differences 
between  these  new  materials. 

AFM  images  recorded  on  freshly  cleaved  NbSe2  crystal  surfaces  exhibit  defects 
typically  7- 10  nm  wide  at  the  surface  and  one  layer  deep  (0.65  ±  0.03  nm).  The  density  of 
defects  is  about  7  x  W'Vnm^.  In  contrast,  freshly  cleaved  M0S2  surfaces  do  not  exhibit 
similar  defects.  Surprisingly,  repetitive  scanning  of  NbSe2  and  M0S2  surface  regions  that 
contain  defects  has  been  shown  to  lead  to  dramatic  wear  and  material  removal.  There  are, 
however,  several  important  differences  in  the  details  of  wear  on  M0S2  versus  NbSe2  (Fig.  3- 
1)  22.29 
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FIGURE  3-1.  Time  sequence  of  images  recorded  on  NbSe2  (upper  panels)  and  M0S2  (lower 
panels). 

First,  for  similar  imaging  forces  (10'^  -  10*^  N)  we  have  shown  that  it  takes  at  least  five  times 
longer  to  wear  a  similar  size  area  on  M0S2  versus  NbSe2.“^ 

Secondly,  we  have  found  that  the  initial  defect  on  the  M0S2  surface  grows 
preferentially  in  a  triangular  shape  from  a  size  of  10  to  >  200  nm  on  edge.  Atomic  resolution 
images  demonstrate  that  wear  occurs  preferentially  along  the  three  equivalent  crystal  lattice 
directions,  and  not  along  the  scan  direction  as  observed  for  NbSe2.  T^he  observations  of 
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highly  selective  wear  and  a  significantly  lower  wear  rate  for  M0S2  indicate  that  this  surface  is 
intrinsically  more  stable  than  that  of  NbSe2.  Notably,  previous  studies  addressing  the 
tribological  properties  of  M0S2  an  NbSe2  did  not  consider  the  higher  stability  of  M0S2  to 
explain  its  superior  propenies.  Our  new  studies  have  shown  that  the  intrinsic  stability  of  the 
M0S2  interface  must  play  an  important  role  in  the  effectiveness  of  this  lubricant. 

3.2.  OXIDATION  OF  M0S2  and  NbSe2 

Since  many  technological  applications  of  transition  metal  dichalcogenide  lubricants 
require  their  use  in  oxidizing  environments  we  have  also  characterized  the  microscopic 
details  of  oxidation  for  M0S2  and  NbSe2. 

AFM  measurements  have  provided  unique  insight  into  this  problem,  and  furthermore, 
these  data  have  also  highlighted  the  distinct  differences  between  M0S2  and  NbSe2.  Oxidation 
of  M0S2  at  440°C  was  found  to  produce  small  defects  in  the  M0S2  surface  as  shown  in 
Figure  3-2. 


FIGURE  3-2.  660  x  660  nm^.  AFM  images  of  (a)  M0S2  and  (b)  oxidized  M0S2.  The  crystal 
was  oxidized  in  air  at  440°C  for  five  minutes. 
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The  size  and  density  of  defects  on  M0S2  following  this  oxidation  are  similar  to  that  found  on 
NbSe2  at  room-temperature.  These  results  also  indicate  that  M0S2  is  intrinsically  more  stable 
than  NbSe2.  Dramatic  support  of  our  hypothesis  has  come  from  oxidation  studies  of  NbSe2. 
Specifically,  oxidation  of  NbScT  at  220°C  for  five  minutes  in  air  causes  extensive  surface 
degradation  as  shown  in  Figure  3-3. 


FIGURE  3-3.  660  x  660  nm^  AFM  image  of  NbSe2  following  oxidation  in  air  at  220°C  for 
five  minutes. 

The  oxidation  of  NbSe2  at  temperatures  «  200°C  less  than  used  for  M0S2  causes  oxide  of 
large  areas  of  the  surface  and  deep  defects  and  contrasts  the  small  defects  observed  in  the 
case  of  MoSo.  These  results  have  thus  clearly  shown  that  M0S2  surfaces  are  more  stable  than 
NbSe2  surfaces.  Furthermore,  since  wear  occurs  preferentially  at  defects,  the  results  from  this 
project  show  that  NbSe2  would  be  a  poor  choice  for  a  lubricant  at  elevated  temperatures  but 
M0S2  may  be  acceptable. 

The  oxidation  of  M0S2  was  also  investigated  over  a  range  of  temperatures  to 
investigate  the  stability  range  of  this  interface  and  the  products  formed  by  oxidation.  AFM 
studies  of  M0S2  oxidized  at  varying  temperatures  for  fixed  time  periods  have  revealed 


several  new  and  important  points.  First,  for  oxidation  temperatures  <  420°C  no  detectable 
defects  are  observed  at  the  M0S2  interface  (Fig.  3-4). 


FIGURE  3-4.  2x2  lam^  AFM  images  recorded  on  M0S2  samples  oxidized  in  O2  for  five 
minutes  at  (A)  400°C  and  (B)  420°C. 

Samples  oxidized  at  >  420°C  show,  however,  defects  and  other  surface  features. 
Oxidation  at  440°C  for  five  minutes  leads  to  the  formation  of  trigonal  pits  20  ±  5  nm  on  edge 
and  one  S-Mo-S  layer  deep,  tr.e  average  density  of  these  defects  is  1.1  x  lO'^/nm^.  As  the 
oxidation  temperature  was  increased  we  found  that  the  defect  density  increased  only  slightly, 
however,  the  size  of  the  defects  increased  significantly  as  shown  in  Figure  3-5. 


FIGURE  3-5.  2  lim  x  2|im  AFM  images  recorded  on  M0S2  samples  that  were  oxidized  in 
O2  for  five  minutes  at  (C)  440°C,  (D)  460°C  and  (E)  480°C. 


After  oxidation  at  460°C  the  triangular  defects  had  increased  in  size  to  approximately  100  nm 
on  edge  and  after  the  480°C  oxidation  these  defects  funher  increased  in  size  to  210  ±  30  nm. 


20 


In  addition,  we  have  also  discovered  the  growth  of  cr>'stallites  on  the  surface  of  M0S2 
following  oxidation  at  temperatures  >  480°C.  The  size  of  the  crystallites,  which  appear  as 
light  plates  in  images,  increases  with  increasing  oxidation  time  and/or  temperature  as  shown 
in  Figure  3-6. 


FIGURE  3-6.  3  |im  x  3  |im  AFM  images  recorded  on  M0S2  samples  that  were  oxidized  at 


480  °C  for  (A)  2,  (B)  5.  (C)  8  and  (D)  1 1  minutes. 
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Analysis  of  the  crystal  surfaces  (following  oxidation)  by  X-ray  photoelectron  spectroscopy 
(XPS)  has  shown  that  the  new  features  detected  in  the  AFM  images  correspond  to  the 
Mo(VI)  oxide,  M0O3.  Hence,  we  have  shown  that  oxidation  of  M0S2  initially  proceeds  (T  > 
440°C)  with  the  formation  of  small  defects,  which  probably  function  as  nucleation  sites,  and 
then  growth  of  micron  size  crystallites  of  M0O3.  This  process  is  shown  schematically  in 
Figure  3-7. 


FIGURE  3-7,  Schematic  view  of  M0S2  before  and  after  oxidation,  [O],  The  M0O3 
crystallite  grows  with  its  b-axis  perpendicular  to  the  MoS2  surface. 

3.3  WEAR  AND  NANOFABRICATION 

Since  it  is  reasonable  to  expect  that  M0O3  also  form  an  M0S2  in  applications 
environments  we  also  characterized  its  tribological  properties  with  the  AFM.  In  general,  we 
have  found  that  both  the  M0O3  and  M0S2  surfaces  are  stable  to  repetitive  scanning  when  the 
imaging  force  is  <  10'^  N.  When  the  imaging  force  is  increased  to  >  5  x  10'^,  however,  the 
M0O3  surface  (but  not  M0S2)  exhibits  wear.  We  have  exploited  this  wear  process  as  a  means 
to  nano-machine  M0O3  in  a  controlled  manner  with  high  resolution.73.27 ,30 
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For  example,  a  line  150  nm  long  is  shown  in  Figure  3-8.  The  width  of  the  lines  in  this 
figure  are  only  10  nm  at  the  M0O3  and  decrease  to  5  nm  at  the  M0O3/M0S2  interface.  The 
average  aspect  ratio  is  thus  on  the  order  of  20:1,  and  it  has  been  possible  to  create  structures 
with  aspect  ratios  as  high  as  40:1.  Such  high-resolution  line  structures  may  have  application 
as  diffraction  gratings. 

In  addition  to  simple  linear  features,  we  have  also  used  the  nanomachining  process  to 
create  significantly  more  complex  structures  (Figure  3-8). 


FIGURE  3-8.  Series  of  500  x  500  nm^  AFM  images  showing  the  machining  of  the  pattern 
"HU".  The  images  were  recorded  nondestructively  using  an  applied  load  of 
1  X  10'^  N;  nanomachining  was  carried  out  using  forces  >  5  x  10'^  N. 

We  have  illustrated  this  point  by  creating  the  pattern  "HU".23.27  ^  important  to  note  that 
the  nanomachining  resolution  does  not  degrade  during  these  operations,  and  thus  it  will  be 
possible  to  use  or  follow  procedures  designed  to  create  structures  with  a  resolution  which 
exceeds  that  obtainable  by  conventional  lithography. 
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We  have  also  used  AFM  for  the  controlled  manipulation  of  nanostructures. 

The  underlying  basis  for  structure  manipulation  in  the  M0O3/M0S2  system  is  that  M0O3  is 
not  strongly  bound  to  the  M0S2  substrate.  It  is  thus  possible  to  pattern  a  structure  in  M0O3, 
separate  this  structure  from  the  M0S2  substrate,  and  then  manipulate  the  object  on  the  M0S2 
surface.  This  process  is  illustrated  in  Figure  3*9  for  a  triangular  object  approximately  80  nm 
on  edge. 


FIGURE  3-9.  Series  of  AFM  images  exhibiting  the  manipulation  of  a  triangular  M0O3 
object  on  the  M0S2  surface.  The  two  translations  shown  in  (b)  and  (c)  each 
correspond  to  a  distance  of  =100  nm. 

It  is  possible  to  move  objects  using  large  forces  (=10'^  N)  and  image  the  resulting 
manipulations  without  further  perturbation  using  a  small  force.  Hence,  we  have  demonstrated 
that  it  is  possible  to  monitor  manipulation  processes  in-situ,  thus  enabling  exquisite  control 
and  precise  positioning  of  objects  and  structures. 
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Section -4 


METAL-SUBSTITUTED  TANTALUM  DISULFIDE 

4. 1  INTRODUCTION 

A  general  and  fundamental  goal  of  our  studies  has  been  to  understand  from  a 
microscopic  perspective  how  dopants  affect  the  properties  of  solids.  To  this  end  we  have  also 
investigated  metal-substitution  in  other  solids,  and  in  panicular.  tantalum  disulfide.^‘-f*  The 
emphasis  in  this  work  has  been  to  probe  how  these  metal-dopants  interact  with  the  change 
density  wave  in  this  material. 

Understanding  the  nature  of  the  interaction  between  impurities  and  a  charge  density 
wave  (CDW)  is  essential  to  understanding  the  static  and  dynamic  propenies  of  the  CDW 
state. 3 1-32  In  general,  pinning  can  be  defined  as  either  strong  or  w.;ak  depending  on  the 
competing  energetics  of  the  CDW-impurity  interaction  and  the  CDW  deformation  energy. 3- 
In  strong  pinning  the  impurity  potential  dominates  me  CDW  elastic  energy  and  pins  the 
phase  of  the  CDW  at  each  impurity  site.  In  weak  pinning  the  CDW  breaks  up  into  constant 
phase  regions  that  are  pinned  collectively  by  impurities.  Despite  the  fundamentally  different 
structural  manifestations  of  strong  and  we; '  oinning,  there  has  been  considerable 
controversy  concerning  the  nature  of  pinning  ..i  CDW  systems. 30-24.33-35  jhis  controversy 
has  been  due  in  large  part  to  the  lack  of  direct  data  characterizing  the  evolution  of  CDW 
structure  with  impurity  doping.  Real-space  imaging  of  the  CDW  structure  in  doped  materials 
by  scanning  tunneling  microscopy  (STM)  is,  however,  beginning  to  provide  a  direct  method 
of  addressing  the  essential  issue  of  strong  versus  weak  pinning.20-24 

In  addition,  an  important  and  general  consequence  of  pinning  is  structural  disorder.  In 
two-dimensional  (2D)  systems  disorder  can  manifest  itself  in  intriguing  ways.  For  example, 
Halperin  and  Nelson  predicted  that  a  2D  solid  can  melt  in  a  continuous  transition  through  a 
hexatic  state  that  is  characterized  by  long-range  orientational  order  and  exponentially 
decaying  positional  order. 36  This  unique  hexatic  state  arises  from  the  formation  of 
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topological  defects,  dislocations,  in  the  lattice.  Disorder  due  to  inipuntv  pinning  has  many 
similarities  to  melting,  although  the  melting  theory  is  based  upon  equilibrium  thermal 
disorder  whereas  pinning  is  typically  a  quenched  disorder.  Since  statistical  averaging  dilfers 
for  equilibrium  and  quenched  disorder-^^  it  is  important  to  examine  carefully  the  analogy  of 
disorder  in  pinned  systems  to  equilibrium  theory. 

We  have  addressed  the  nature  of  CDW  impurity  pinning  in  Nb-doped  TaSs, 
NbxTai.xS2,  and  have  elucidated  the  detailed  characteristics  of  disorder  in  this  pinned  CDW 
system.  Real-space  images  of  the  CDW  lattice  as  a  function  of  iNb  impurity  concentration 
were  obtained  by  STM,  and  were  quantitatively  analyzed  to  determine  the  topological  defects 
and  the  radial,  translational  and  orientational  correlation  functions.  These  data  show 
unambiguously  that  CDW  pinning  by  the  Nb  impurities  is  weak  or  collective.  In  addition,  our 
work  demonstrates  that  this  system  evolves  through  crystalline,  hexatic  glass  and  amorphous 
.states  as  a  function  of  impurity  concentration. 

4.2  STM  CHARACTERIZATION 

High  resolution  STM  images  of  the  NbxTai-xS2  materials  are  shown  in  Figures  4-1  to 
4-3.  Images  of  the  incommensurate  CDW  phase  of  undoped  TaS2  samples  exhibit  a  well- 
ordered  hexagonal  CDW  superlattice  and  atomic  lattice  (Fig.  4-1). 
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FIGURE  4- 1 .  STM  image  of  a  TaS2  recorded  at  360  K.  The  image  exhibits  both  the  atomic 
lattice  (a  =  3.35  A)  and  CDW  super  lattice  (a  =  1 1.8  A). 

In  contrast,  substitution  of  Nb  causes  disorder  in  the  CDW  lattice  (Figs.  4-2,  4-3).  The 
images  of  the  x(Nb)  =  0.02  samples  exhibit  areas  in  which  the  CDW  lattice  has  hexagonal 
order,  and  also  regions  containing  defects;  these  defects  introduce  disorder  into  the  CDW 
lattice. 
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FIGURE  4-2.  STM  images  of  the  incommensurate  CDW  state  in  (a)  Nb0.02Ta0.98S2  and 
(b)  Nb0.04Ta0.96S2.  Black  lines  highlight  the  insertion  of  an  extra  row  of 
CDW  lattice  sites  in  (a).  Two  distinct  loops  and  Burgers  vectors  are  drawn  in 
(b)  to  highlight  dislocations  in  the  CDW  lattice.  A  5-fold-7-fold  disclination 
pair  is  also  shown  in  (b). 


The  predominant  defects  observed  in  the  samples  containing  x(Nb)  <  0.04  are  dislocations. 
Dislocations  are  formed  by  the  insertion  of  an  extra  half  row  of  CDW  sites  in  the  lattice; 
black  lines  in  Figure  4-2  highlight  the  creation  of  one  dislocation.  An  important  point  to 
recognize  about  these  defects  is  that  there  is  a  significant  strain  field  at  the  dislocation. 38  The 
CDW  can  relax  the  strain  field  by  locally  deforming  (i.e.,  rotating),  although  such  local 
rotations  create  disorder. 

As  the  impurity  concentration  increases  to  x(Nb)  =  0.04  we  find  that  dislocations 
appear  with  a  higher  density  than  in  the  x(Nb)  =  0.02  samples  (Fig.  4-2b).  The  CDW  lattice 
rows  near  the  dislocation  are  deformed  as  discussed  above.  However,  in  areas  free  of 
dislocations  the  CDW  lattice  is  locally  ordered.  We  have  highlighted  several  dislocation 
cores  (the  region  where  the  CDW  deforms)  by  constructing  Burgers  loops  (Fig.  4-2b).  The 
Burgers  loop  consists  of  an  equal  number  of  steps  along  each  lattice  direction;  the  loop  will 
remain  open  if  it  encloses  a  single  dislocation.  The  vector  pointing  from  start  to  end  of  the 
loop,  the  Burgers  vector,  uniquely  defines  the  dislocation. 38  We  find  that  the  Burgers  vectors 
defining  the  CDW  dislocations  in  the  Nb-doped  samples  occur  along  each  of  the  three 
crystallographic  axes,  and  thus  we  conclude  that  impurity-induced  dislocations  occur 
randomly  in  the  CDW  lattice.  In  the  samples  containing  higher  impurity  concentrations, 
x(Nb)  =  0.07  and  0. 10  the  STM  images  exhibit  extended  defects  (Fig.  4-3). 
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FIGURE  4-3.  STM  images  of  the  incommensurate  CDW  phase  in  (a)  Nbo.07Tao.93S 2  and 
(b)  Nbo.1Tao.9S2. 

The  CDW  lattice  in  these  latter  samples  exhibit  significant  disorder  with  regions  of 
hexagonal  order  extending  only  several  lattice  constants. 

To  elucidate  in  detail  the  topology  and  density  of  these  defects  we  have  quantitatively 
analyzed  the  STM  images.  The  quantitative  analysis  involves  defining  the  x,y  coordinates  of 
each  CDW  maxima  and  the  unique  nearest  neighbors  in  the  lattice.  The  CDW  maxima. 
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which  we  call  lattice  sites,  were  located  to  ±1  pixel  accuracy  by  a  thresholding  algorithm. 
Once  the  lattice  points  are  located,  a  sweepline  algorithm  is  used  to  construct  the  Voronoi 
diagram  for  the  lattice.  The  Voronoi  diagram  determines  the  nearest  neighbors  of  all  of  the 
lattice  points  uniquely, and  thus  can  be  used  to  illustrate  explicitly  all  topological  defects  in 
the  lattice.  To  illustrate  these  defects  we  triangulate  the  Voronoi  diagram  by  drawing  "bonds  " 
from  all  CDW  lattice  points  to  their  nearest  neighbors.  Hence,  fully  coordinated  lattice  sites 
are  indicated  by  zIk  bonds  to  the  CDW  maxima,  while  defects  contain  fewer  or  greater  bonds. 
Typical  results  obtained  from  the  analysis  of  x(Nb)  =  0,  0.02,  0.04,  0.07  and  0. 10  samples  are 
shown  in  Figure  4-4. 
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FIGURE  4-4.  Triangulations  of  the  STM  images  recorded  on  NbxTai.xS2  crystals  where  (a) 
.V  =  0,  (b)  .r  =  0.02,  (c)  x  =  0.04,  (d)  x  =  0.07,  and  (e)  x  =  0.10.  Lattice  sites 
that  do  not  have  six-fold  coordination  are  highlighted  by  shading  the 
trianeles  which  have  the  non-six  fold  coordinate  site  as  a  vertex. 


In  these  triangulations  we  have  highlighted  the  defect  (non-six-fold  coordinate)  sites  by 
shading.  Analy.ses  of  images  recorded  on  pure  TaS2  show  that  the  CDW  lattice  is  free  of 
topological  defects;  that  is,  all  of  the  lattice  sites  are  sixfold  coordinate  (Fig.  4-4a).  The 
triangulations  explicitly  show,  however,  that  the  Nb-doped  materials  have  topological  defects 
in  the  CDW  lattice.  At  low  concentrations  of  impurities,  x(Nb)  =  0.02  and  0.04,  we  find  that 
the  dislocations  consist  of  fivefold/sevenfold  disclination  pairs  (Fig.  4-4).  Extended  defect 


networks  are  also  obvious  in  the  triangulation  data  for  the  x(Nb)  =  0.07  and  0. 10  samples. 
These  extended  topological  defects  consist  of  dislocations  and  free  disclinations. 

We  have  also  used  this  data  to  determine  the  separation  between  dislocations.  We  find 
that  the  average  spacing  between  dislocations  in  the  x(Nb)  =0.02, 0.04,  0.07  and  0. 1  CDW 
lattices  is  12,  8,  5,  3  lattice  constants,  respectively.  In  comparison,  the  average  separation 
between  the  Nb  impurities  in  the  lattice  (dj)  is  0.80, 0.57, 0.43  and  0.36  CDW  lattice 
constants,  respectively.  The  average  separation  between  dislocations  is  thu^  cd^'^ays  greater 
than  the  average  impurity  spacing;  that  is,  impurity  pinning  of  the  CDW  must  be  a  collective 
effect.  Hence,  these  results  show  that  CDW  pinning  in  the  NbxTai.xS2  materials  is  weak. 

4.3  THEORETICAL  ANALYSIS 

To  further  examine  the  issue  of  CDW  pinning  we  consider  theoretical  models  of  the 
CDW-impurity  interaction.  The  one-dimensional  model  first  proposed  by  Fukuyama,  Lee 
and  Rice  (FLR)  has  been  the  most  widely  studied  example  for  CDW-impurity  pinning.32  For 
this  one-dimensional  model  the  CDW  wave  function  is 

p(x)  =  po  -t-  pcos(Q-x  +  (l)(x))  ( 1 ) 

where  p(x)  is  modulated  charge  density,  po  is  the  background  charge  density,  p  is  the  CDW 
amplitude,  Q  is  the  CDW  wave  vector,  and  (j)(x)  is  CDW  phase.  Analysis  of  this  model  shows 
that  there  are  two  distinct  regimes  which  describe  the  strength  of  the  CDW-impurity 
interaction  (i.e.,  strong  and  weak  pinning).  These  two  regimes  can  be  defined  in  terms  of  a 
competition  between  the  pinning  energy  Ep  =  VoP  and  the  CDW  elastic  energy  Eg  =  aK.  In 
strong  pinning  the  impurity  potential  Vq  dominates  the  elastic  energy:  e  =  Vop/ttK  »  1. 

Since  the  CDW  is  pinned  to  each  impurity  site  in  the  lattice  the  CDW  phase  undergoes 
abrupt  changes  at  the  impurity  sites.  Conversely,  in  weak  pinning  the  CDW  stiffness  resists 
deformation  by  the  impurity  potential:  e  «  1.  For  weak  pinning  the  total  energy  is 
minimized  through  collective  interactions.  Hence,  the  CDW  phase  can  vary  smoothly  over  a 
region  containing  many  impurities. 
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The  2D  CDW  system  studied  here  should  form  areas  with  local  phase  coherence 
and  positional  order  if  the  system  is  weakly  pinned.  Dislocations  and  other  defects  will 
separate  phase  coherent  regions.  These  locally  coherent  regions  can  be  understood  by  scaling 
arguments.  For  an  area  of  length  scale  L  with  impurity  density  nj,  the  energy  gain  from 
pinning  per  unit  area  will  be  Vop(L2ni)*/2L-2  while  the  elastic  energy  cost  is  -aicL  -.'^O  By 
minimizing  the  total  energy  density  one  finds  an  optimized  correlation  length  Lc, 

Lc  =  (2/e)di  (2) 

Where  di=(l/ni)l/2  is  the  spacing  between  impurities.  Within  the  length  scale  Lc,  the  CDW 
phase  varies  smoothly  between  -n  and  K.;  positional  correlations  will  only  be  sustained  on  a 
scale  smaller  than  Lc-  Alternatively,  this  scaling  argument  can  be  used  to  explain  the 
formation  of  dislocations.  The  energy  cost  for  a  dislocation  loop  with  size  L  is  ~aic  up  to  a 
logarithm,38  while  the  energy  gain  due  to  pinning  is  VopCL^ni)!/^.  Hence,  creating  a 
dislocation  on  length  scale  larger  than  Lc=(2/e)di  is  energetically  favorable. 

As  discussed  above  the  spacing  between  dislocations  in  the  CDW  lattice  is 
significantly  larger  than  the  distance  between  Nb  impurities.2'^  This  observation  is  consistent 
with  the  weak  pinning  model  developed  above.  A  related  criteria  for  weak-pinning  is  the 
presence  of  ordered  regions  with  a  critical  size  directly  proportional  to  the  impurity  spacing. 
Notably,  analysis  of  large  area  STM  images  demonstrate  that  areas  of  the  CDW  lattice  that 
are  free  from  dislocations  exhibit  a  smooth  variation  in  the  CDW  phase  as  we  'xpect  for  the 
case  of  weak  pinning.  Moreover,  we  can  use  our  scaling  equation  to  estimate  the  magnitude 
of  the  dimensionless  parameter  e.  Using  the  Lc  determined  from  our  images  and  the  known 
impurity  spacing  we  find  that  e  =0. 1-0.2  «  1.  These  results  demonstrate  convincingly  that 
the  CDW  pinning  in  NbxTai.xS2  is  weak. 

It  is  also  important  to  examine  the  consequences  of  the  disorder  produce  by  weak 
pinning  in  this  system.  We  first  examine  the  structure  factor,  S(k)  =  lp(k)l2,  to  investigate  the 
average  structural  effects  of  pinning.  S(k)’s  were  calculated  from  the  square  of  the  Fourier 


transform  of  the  STM  images;  the  results  for  the  x{Nb)  =  0.02,  0.04,  0.07  and  0. 10  samples 
are  shown  in  Figure  4-5. 


FIGURE  4.5  Structure  factors  corresponding  to  the  CDW  lattices  of  (a)  x  =  0.02,  (b)  x  ~ 
0.04,  (c)  X  =  0.07,  and  (d)x  =  0.10  NbxTai.xS2  materials. 

S(k)  for  pure  TaS2  (not  shown)  exhibits  sharp  sixfold  symmetric  peaks.  These  peaks  broaden 
both  radially  and  angularly  as  the  impurity  concentration  increases  to  0.04.  For  the  x(Nb)  = 
0.07  and  0.10  the  first  order  Bragg  peaks  have  broadened  to  form  a  ring  whose  intensity  has  a 
sixfold  modulation.  The  angular  broadening  is  due  to  CDW  rotations  that  were  discussed 
above  in  section  A.  This  broadening  indicates  a  loss  of  orientational  order;  however,  it  is  not 
possible  to  provide  a  quantitative  measure  of  the  orientational  disorder  from  this  S(k)  data.  In 
addition,  the  radial  widths  of  the  Bragg  peaks  indicate  that  the  translational  correlation  length 
is  short  in  all  of  the  Nb-doped  materials.  Notably,  these  S(k)  data  are  qualitatively  similar  to 
the  results  expected  for  2D  melting  through  a  hexatic  state.^^ 
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To  provide  a  quantitative  measure  of  the  order  in  this  system  and  how  it  is  affected  by 
impurity  pinning  we  have  investigated  the  radial  distribution,  translational  correlation  and 
orientational  correlation  functions.  The  radial  distribution  function  is  defined  as  g(r)  = 
<n(r)>/No,  where  n(r)  is  the  point  density  at  distance  r  from  the  origin  of  the  structure  and  Nq 
is  the  average  density  of  points.  For  a  general  structure,  g(r)  tends  to  1  as  r  ->«>;  that  is,  the 
distribution  of  points  at  large  distances  appears  uniform.  The  results  of  the  calculations  of 
g(r)  are  presented  in  Figure  4-6. 
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x  =  0.0 


r 

FIGURE  4-6.  Radial  distribution  functions  for  the  .r  =  0,  x  =  0.02,  x  =  0.04,  x  =  0.07,  and  .v 
=  0. 10  NbxTai.xS2  materials.  The  x  axis  corresponds  to  CDW  lattice 
constants. 

For  the  x(Nb)  =  0  sample  g(r)  oscillates  strongly  out  to  more  than  10  lattice  constants.  The 
persistent  oscillation  in  g(r)  indicates  that  positional  order  is  long  range.  When  the  impurity 
level  increases  to  0.02  and  0.04  the  amplitude  of  the  oscillations  in  g(r)  decreases  and  the 
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oscillations  die  out  much  more  rapidly  than  in  g(r)  for  the  pure  material.  For  the  x(Nb)  =  0.07 
or  0. 1  materials  the  amplitude  of  the  peaks  in  g(r)  is  funher  reduced  and  these  peaks  only 
persist  to  the  second  or  third  nearest  neighbors.  These  results  indicate  that  Nb  impurities 
destroy  translational  order  in  this  CDW  system. 

To  further  characterize  the  order  of  this  2D  system  we  have  calculated  the 
translational  and  orientational  correlation  functions.  Previous  theoretical  work  has  shown  that 
these  functions  provide  a  unique  description  of  the  order  in  2D  systems.^^  The  translational 
correlation  function  is  defined  as  GT(r)  =  where  \|r(r)  is  the  translational 

3 

order  parameter  at  position  r:  y(f)=  ^  exp(iGi  •  r) /  3.  The  summation  is  over  the  three 

1=1 

reciprocal  lattice  vectors,  Gl-  The  orientational  correlation  function  is  G(5(f )  = 
■^V6(0)Vd(r)*^'  where  vi/6(r)  ‘S  the  orientational  order  parameter: 

Ni 

\|/6=  ^  exp(i60(r))  /  Ni .  6(1)  is  the  bond  angle  with  respect  to  a  fixed  direction,  where  the 

lal 

"bonds"  are  the  lines  that  connect  nearest  neighbors  as  in  the  triangulation  (Fig.  4).  Nj  is  the 
number  of  bonds  at  the  i‘^  point  and  the  summation  is  over  all  bonds  at  position  r.  For  a 
perfect  lattice  GT(r)  =  G6(f)  =  1. 

Dislocations  observed  in  the  CDW  lattice  will  reduce  both  GT(r)  and  G6(r), 
although  the  detailed  effect  on  these  functions  is  expected  to  be  quite  different.  Qualitatively, 
a  dislocation  is  veiy  effective  in  destroying  translational  order  since  it  displaces  the  CDW 
sites  on  the  order  of  one  lattice  constant.  Translational  order  should,  therefore,  only  be 
sustained  in  regions  free  of  dislocations.  Indeed,  Halperin  and  Nelson  showed  that  when  the 
structural  order  is  dominated  by  dislocations  the  translational  correlation  length  *  the 
dislocation  spacing.  In  the  case  of  G6(r)  the  effect  of  dislocations  is  expected  to  be  quite 
different.  For  a  dislocation  consisting  of  a  bound  pair  of  fivefold-sevenfold  disclinations  (Fig. 
4),  the  bond  orientational  order  parameters  at  these  two  points  are  obviously  reduced. 

Because  the  disclinations  are  closely  bound,  they  do  not  change  the  coordination  number  and 
orientation  of  points  away  from  the  dislocation.  Hence,  dislocations  do  not  necessarily 
destroy  orientational  order. 
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The  calculated  results  for  GtI? )  and  G6(r)  are  shown  in  Figure  4-7. 
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FIGURE  4-7.  Translational  correlation  functions  Gjir)  (a-e)  and  orientational  correlation 
functions  GgCr)  (f-j)  calculated  from  images  x  =  0,  0.02, 0.04,  0.07,  and  0. 10 
NbxTai.xS2  materials.  The.r  axes  in  these  figures  correspond  to  lattice 
constants.  The  points  represent  the  experimental  data,  and  the  lines  correspond 
to  exponential  fits  (a-e)  and  (i,  j)  or  power  law  (r-H)  fits  (g,  h). 

GT(r)  and  G6(r)  decay  very  little  over  20  lattice  constants  for  the  pure  sample.  However,  all 
of  the  samples  that  contain  Nb  impurities  exhibit  a  rapid  decay  of  translational  order.  This 
decay  in  GT(f )  can  be  fit  reasonably  well  to  an  exponential  of  the  form  GT(r)  =  expC-^j/r). 
where  is  the  translational  correlation  length.  The  range  of  the  ^t's  determined  by 
averaging  a  number  of  images  are  ^7  =  7-10,  3-6.  2-3,  and  1-2  lattice  constants  for  the  x(Nb) 
=  0.02,  0.04,  0.04,  and  0. 10  samples,  respectively.  The  values  of  ^7  are  similar  to  the  average 
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dislocation  spacing  determined  from  triangulation  data  (e.g..  Fig.  4-2),  and  funhermore  these 
results  demonstrate  quantitatively  that  weak  pinning  by  the  Nb  impumies  destroys 
translational  order.  In  contrast,  we  find  that  the  orientational  correlations  die  out  slowly  for 
the  0.02  and  0.04  samples,  although  G^Cf )  decays  rapidly  for  x(Nb)  =  0.07  an  i  0. 10  (Figs. 
8f-j).  If  the  x(Nb)  =  0.02,  0.04,  0.07  and  0.10  data  are  fit  with  an  exponent!,.  C(,(r)  =  exp(- 
^6/r),  we  obtain  correlation  lengths  of  ^6  =  200,  100,  1 1,  and  5  lattice  constants,  respectively. 
A  better  fit  to  G^fr)  for  x(Nb)  =  0.02, 0.04  is  obtained,  however,  using  an  algebraic  decay, 
G6(f )  =  r-h,  with  T)  =  0.03  and  0. 12,  respectively. 

The  simultaneous  observation  of  long  range  orientational  order  and  short  range 
translational  order  in  the  x(Nb)  =  0.02  and  0.04  samples  is  strongly  suggestive  of  the  hexatic 
state  that  arises  in  2D  melting  driven  by  topological  defects.36  2D  melting  theory,  vvhicn  is 
based  on  equilibrium  statistical  mechanics,  predicts  that  melting  is  a  continuous  phase 
transition  between  (1)  a  crystalline  solid  pha.se  that  has  long  range  translational  and 
orientational  order:  Gj((F)  ~  r'^,  G^(r)  ~  1;  (2)  a  hexatic  phase  that  has  short  range 
translational  order  but  long  range  orientation  order:  Gj(r)  =  e''’^T,  G^Y)  -  r’’^;  and  (3)  a 
liquid  phase  that  has  short  range  translational  and  orientational  order:  Gj(f )  =  e''^^^T,  G^r)  = 
e'^'^^6.  Our  data  show  similar  behavior  and  thus  we  suggest  that  the  CDW  lattice  evolves  from 
crystalline  (x=0)  to  hexatic  (x=0.02-0.04)  to  amorphou*  ,x>0.07)  stales. 
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Section-5 


Summary 

Scanning  tunneling  microscopy  (STM)  and  atomic  force  microscopy  (AFM)  have 
been  used  to  characterize  the  atomic  level  structure  of  electronic  properties,  reactivity  and 
wear  of  metal  dichalcogenide  materials  that  are  or  have  potential  as  solid  state  lubricants. 
Single  crystals  of  M0S2,  NixMoi.xS2,  MoS2-xSex  and  MoS2-xTex  have  been  prepared  to 
determine  how  chemical  modifications  affect  the  local  structure  and  electronic  properties  of 
this  lubricant.  STM  images  of  Ni-doped  M0S2  show  localized  electronic  states  due  to  the  Ni 
atoms,  while  images  of  Se-  and  Te-doped  materials  indicate  that  anion  substitution  is 
electronically  delocalized.  AFM  studies  of  Te-doped  M0S2  show,  however,  that  the  tellurium 
doptmts  form  atomic  scale  structural  protrusions  that  may  reduce  sliding  friction. 

AFM  has  also  been  used  to  characterize  nanometer  scale  wear  and  oxidation  on  M0S2 
and  NbSe2  surfaces.  In  atmosphere  at  room- temperature  AFM  studies  showed  that  NbSe2 
wears  approximately  three  times  faster  than  M0S2.  Funhermore,  oxidation  studies 
demonstrated  that  NbSe2  was  significantly  more  reactive  than  M0S2  with  molecular  oxygen. 
These  results  indicate  that  the  intrinsic  stability  of  the  M0S2  surface  make  it  an  effective 
lubricant.  AFM  was  also  used  to  elucidate  the  growth  of  M0O3  on  the  surface  of  M0S2 
during  oxidation,  and  to  study  wear  properties  of  these  M0O3  crystallites.  The  AFM  tip  was 
used  to  define  lines  with  10  nm  resolution  in  M0O3  and  to  manipulate  distinct  M0O3 
structures  on  the  M0S2  surface. 

In  addition,  metal-substitution  in  TaS2  has  been  studied  systematically  using  STM 
and  theoretical  methods.  Investigations  of  Nb-doped  materials  provide  the  first  direct 
structural  evidence  for  weak  pinning  of  an  electronic  lattice,  and  furthermore,  theoretical 
analyses  have  led  to  the  discovery  of  a  new  phase,  the  hexatic,  in  these  materials. 
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